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We study the stability and synchronization of predator-prey populations subjected to noise. The 
system is described by patches of local populations coupled by migration and predation over a 
neighborhood. When a single patch is considered, random perturbations tend to destabilize the 
populations, leading to extinction. If the number of patches is small, stabilization in the presence 
of noise is maintained at the expense of synchronization. As the number of patches increases, both 
the stability and the synchrony among patches increase. However, a residual asynchrony, large 
compared with the noise amplitude, seems to persist even in the limit of infinite number of patches. 
Therefore, the mechanism of stabilization by asynchrony recently proposed by R. Abta et. al.0], 
combining noise, diffusion and nonhnearities, seems to be more general than first proposed. 



PACS numbers: 87.23.Cc,05.45.Xt,87.18.Hf 

The model proposed independently by Lotka [l[ and 
Volterra was probably the first to describe mathemat- 
ically the dynamics of predators and preys. Its success 
and widespread use by early biologists is mostly due to 
its ability to qualitatively describe the population oscil- 
lations of both preys and predators @ . The model, how- 
ever, is well known to be unstable under the addition of 
noise, which causes the amplitude of the population oscil- 
lations to increase until one (or both) species eventually 
becomes extinct. Stability, in the sense of co-existence of 
both species, can be regained if several patches of pop- 
ulations are coupled via dispersal or predation. Compu- 
tational simulations have shown that stability increases 
(with respect to noise amplitude) with the number of 
patches considered 0, H, [|[ . The ultimate reason for the 
stabilization was recently pinned down by R. Abta et. 

, who studied in detail the case of two patches. They 
concluded that the crucial condition for stabilization is 
the development of an asynchrony between the popula- 
tion oscillations in each patch, resulting from the com- 
bined action of diffusion and noise. Such an asynchrony 
develops if the frequency of the population oscillations 
depends on their amplitude. 

Although the results of rcf.[3] are of theoretical and 
conceptual importance, simple Lotka- Volterra equations 
(LV, for short) are seldom used to model population dy- 
namics now-a-days. Instead, models involving logistic 
type of interactions displaying limit cycles or chaotic at- 
tractors have become common 

@, M M El- In this 
report we discuss the problem of stabilization and syn- 
chronization for a predator-prey system displaying an at- 
tracting limit cycle Q. For this system all asymptotic 
orbits have the same oscillation frequency, which is the 
frequency of the attractor. As a consequence, orbits dis- 
placed from the attractor by a perturbation will return 
to the attractor and to the same frequency of oscillation, 
different from LV model. We shall call this system LC 
(limit cycle) for short. The system is only weakly sta- 
ble in the presence of noise, in the sense that it becomes 
unstable when the noise amplitude crosses a threshold 
which is very small. Therefore, for noise amplitudes that 



are not too small, LC behaves like LV, with noise driving 
one or both species to extinction. 

A spatial version of the LC model, SLC, can be con- 
structed by allowing patches of local populations to in- 
teract. If the patches are strongly coupled, the dynamics 
in each spatial region synchronizes and the system be- 
haves like a single well mixed population, identical to 
the original LC model. Here we study the stability and 
synchronization of the spatial model under random per- 
turbations in the strongly coupled regime. Because of the 
attracting limit cycle, it is not clear that an asynchrony 
among patches will develop, since perturbed orbits will 
always have similar amplitudes and, therefore, similar os- 
cillation frequencies. For a system with only two patches, 
we show that desynchronization indeed takes place, lead- 
ing to the stabilization of the population oscillations, in 
agreement with the results of Q for the LV model. As 
more patches are added, SLC system becomes stable un- 
der larger noise amplitudes [J] and the asynchrony de- 
creases exponentially with system size. However, a resid- 
ual asynchrony, much larger than the noise amplitude, 
seems to survive even in the limit of infinite number of 
patches. Therefore, even for large systems displaying an 
attracting limit cycle, the combined action of diffusion 
and noise still plays a crucial role in desynchronizing the 
patches. 

The LC predator-prey model with noise is given by the 
equations 
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Vn+i = Vn[e dl + F y (x n )] + rjy 



where P x {y) = e~ y l a accounts for the predation of y 
upon the prey x and F y (x) = 1 — e~ x ^ for the reproduc- 
tion of y due to the feeding upon x. The variables rj x and 
r\ y are random numbers representing the external noise 
whose distributions are homogeneous and limited to the 
interval [—A, A]. 

To understand the role of each term in these equations 
let us first consider r\ x = r\ y = 0. Then, in the absence 
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of predators, P X (Q) = 1 and the population of preys con- 
verges to the normalized value x = 1, provided a < 1. 
As the number of predators increases. P x {y) decreases, 
reducing the population of preys. Similarly, in the ab- 
sence of preys F y (0) = and the population of predators 
decreases steadily because of the intrinsic death rate d\. 
When the number of preys is sufficiently large so that 
g-di _|_ F y (x) > 1 the number of births becomes larger 
than the number of deaths and the population of preda- 
tors grows. 

Following Q we fix the model parameters at a = e _1 , 
a = 0.8 and (3 = 0.4. In the absence of noise the system 
displays two kinds of attractive orbits, depending on the 
value of the death rate d\\ for < d\ < 0.6 the attractor 
is a limit cycle and for d\ > 0.6 it is a fixed point. When 
the noise is turned on it might happen that the popula- 
tion densities x or y become less than zero. In this case 
we set it back to zero to avoid negative values. 

Throughout this paper we have fixed the death rate 
of the predator at d\ = 0.1. In this case the trajecto- 
ries converge to a limit cycle where the prey population 
oscillates between 2.7 x 10~ 4 and 0.3. This attractor is 
only very weakly stable, since the addition of noise with 
amplitudes of the order of A = 6 x 10~ 5 is enough to 
drive both species to extinction. 

To consider more than one patch we extend the LC 
model to M: 
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where i,j label the position of the patches on a two- 
dimensional grid of N x N patches. This spatial version 
(SLC) allows migration of predators and preys, with rates 
m y and m x respectively, and predation over a predation 
neighborhood R. Here we consider square neighborhoods 
with sides 2R + 1 centered on the predator. 

The sum over I and m in the migration terms is 
restricted only to the four nearest neighbors of the 
site The averages (P x (y n )) R are performed 

over the patches that are in predation neighborhood: 
(P^(y)) R = iEL-K^(l/ ,+iJ+m ) where N R = 
(2R + l) 2 is the number of patches within the neighbor- 
hood. The feeding function F y ((x) Ri:j ) is calculated over 
the average number of x on the predation neighborhood: 



F y ((x) R ) = F y (^EL=-^ +U+m j- We fix the mi- 
gration rate of the prey and the predator at m x = 0.01 
and m y = 0.1 respectively, so that the migration rate of 
the predator is larger then that of the prey [HI, [l4| • 



In the SLC model, the average dynamic behavior can 
be very different from that of the LC model. However, 
if the predation neighborhood encompasses most of the 
patches in the grid, the system becomes strongly coupled 
and the time evolution in each patch synchronizes with 
the others 

Two patches 

In our first analysis we consider only two patches. In 
this case Eq.([2|) reduces to 
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where i labels one of the patches and j the other one. 
The averages are given by (P x (y n )) = \{Px{y l ) + Px{y 3 )) 
and Fy((x n )) = F y ( x "t x )■ In the absence of noise and 
with random initial conditions, the time evolution of the 
populations synchronize perfectly. However, when noise 
is added to the system, the populations oscillate between 
synchronized and desynchronized phases. In our simu- 
lations we let the system evolve without noise for 2,000 
time steps, which is enough to synchronize the patches. 
Noise with amplitude of A = 6 x 10~ 5 (which is enough 
to destabilize the dynamics on a single patch) is then 
added to the system, which is further evolved for another 
28,000 steps. Figure Q] shows the populations in one of 
the patches in the phase-space for times between 1,000 
and 10,000. The coupling between the patches stabilizes 
the system. For larger noise amplitudes, of the order of 
A = 2 x 10~ 4 , instability sets in again and both species 
may go extinct. 

In order to quantify the synchrony between the 
patches, we fixed a reference point at r"o = (0.04, 0.75) 
and measured the phase and amplitude differences 
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2 , where fi is the phase space 
position of the populations in patch i. Figure [21 shows 
the phase difference \0qi — #02! and amplitude difference 
|?*oi — tq2 |. As the initial conditions for each patch 
arc different, the phase and amplitude differences are 
initially non-zero. However, the patches quickly syn- 
chronize in the absence of noise (first 2,000 time steps) 
and de-synchronize again when noise is added. 

Increasing the number of patches 

When more patches are taken into account, the sta- 
bility properties change qualitatively. We considered the 
spatial model on a set of N x N patches, Eq. © , with pe- 
riodic boundary conditions. As before, we fixed d\ =0.1, 
a = 0.8, a = 0.4, m x = 0.01, m y = 0.1 and A = 6 x 10~ 5 . 
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Figure 1: Phase-space trajectories for the case of two patches. 
The plots show the time steps between 1,000 and 10,000 for 
only one patch. The patches are synchronized up to 2,000 
steps, when noise is added and synchronization breaks down. 
The patches nearly synchronize when x ~ and get signifi- 
cantly desynchronized for large x. 
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Figure 2: Phase (in radians) and amplitude differences be- 
tween the two patches for the trajectories shown in Figure fT}. 
The crests correspond to regions where x is large. 



and calculated the phase and amplitude differences for 
each patch with respect to patch number one. The 
temporal averages of the phase and amplitude differ- 
ences were computed for each individual patch for the 
last 50,000 time steps, {\9\ — 0i\) and (|ri — fi\) and 

also the global average, A9 = N 2_i 2i=2 (l^i — an< ^ 

Af=^E£2<l*Wil>- 

In agreement with D. D. Donalson et. al. [J we found 
that the populations resist to higher noise amplitudes 
as the number of patches increases. Moreover, for fixed 
noise amplitude and coupling ratio, the patches tend to 
become more synchronized as TV increases. Therefore, 
for large systems, stability and asynchrony are not as 
correlated as in the case of two patches Q . 

Figure [3] shows that the average asynchrony decreases 
exponentially with grid size, following approximately the 
curve y(N) = Aexp(—N/b) + c. For the amplitude dif- 
ference (gray curve) we found A = 0.017, b = 18 and 
c = 0.016 and for the phase difference A = 0.011, b = 16 
and c = 0.0046. In the limit of infinitely many patches 
the average amplitude difference tends to c = 0.046, 
which is significantly larger than the noise, A = 6 x 10 -5 . 
The conclusion of these numerical experiments is that the 
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Figure 3: Global average of phase difference A8 (black), and 
amplitude difference Af (gray), as a function of grid size, TV. 
The lines are first order exponential decay fit. 



The only free parameters are the sizes of the grid, TV, and 
the predation radius, R. As in the case of two patches, 
it is possible to synchronize the patches in the absence of 
noise if the coupling is sufficiently strong Q. 

In order to study the effect of noise on large synchro- 
nized systems we did simulations in which ./V was varied 
but the ratio (2R + 1)/N, was kept constant. We consid- 
ered 8 combinations of TV and R with (2R+1)/N = 5/6: 
(1) TV = 6, R = 2; (2) TV = 18, R = 7; (3) TV = 30, 
R = 12; (4) TV = 42, R = 17; (5) TV = 54, R = 22; (6) 
TV = 66, R = 27; (7) TV = 78, R = 32; (8) TV = 90, 
R = 37. In all cases the patches synchronized in the 
absence of noise. 

The simulations started with random initial conditions 
and were iterated by 60,000 time steps. Noise was added 
only after the first 5,000 steps. As in the case of two 
patches, we fixed the reference point ro = (0.04,0.75) 



simultaneous presence of noise and diffusion seems to lead 
to significant asynchrony even if the unperturbed dynam- 
ics has an attractive limit cycle. For small systems this 
asynchrony can stabilize the populations by allowing the 
migration of individuals from more populated patches to 
those where extinction is imminent. For large systems, 
asynchrony decreases exponentially fast with system size 
but never disappears completely. 

The mechanism responsible for the desynchroniza- 
tion in systems with attractors seems to be the same 
proposed in Q, i.e., motion with amplitude depend 
frequencies. Since the dynamics tend to bring perturbed 
orbits back to the attractor, a typical trajectory always 
wanders in the vicinity of the attractor. However, the 
nonlinear character of the equations amplify these small 
deviations producing significant frequency differences 
that are reflected in the desynchronization. Therefore, 
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the mechanism of R. Abta et. al.[7|, combining noise, 
diffusion and nonlinearities, seems to be more general 
than first proposed. 



Acknowledgements This work was partially supported 
by Fapcsp and CNPq. 



[1] A.J. Lotka Elements of Physical Biology, Baltimore, 

Williams and Willkins (1925). 
[2] V. Volterra, Journal du Conceil, Conceil International 

pour l'exploration de la Mer, 3 3 (1928). 
[3] J. D. Murray. Mathematical Biology. Springer, Berlin, 

1993. 

[4] D. D. Donalson and R. M. Nisbet, Ecology. 43, (1999). 
[5] W. G. Wilson, A. M. Deroos and E. Mccauley, Theor. 

Popul. Biol 43, (1993). 
[6] C. J. Brigges and M. F. Hoopes, Theor. Popul. Biol 65, 

(2004). 

[7] R. Abta, M. Schiffer and N. M. Shnerb, Phys. Rev. Let- 
ters 98, 098104 (2007). 

[8] S. B. L. Araujo and M. A. M. de Aguiar, Phys. Rev. E 
75, 061908 (2007). 



[9] W. M. Schaffer, The American Naturalist. 124, 798 
(1984). 

[10] A. Hastings and T. Powell, Ecology. 72, 896 (1991). 
[11] A. Hastings, C. L. Horn, S. Elmer, P. Turchin and H. C. 

J. Godfray, Annual Review of Ecology and Systematics. 

24, 1 (1993). 

[12] J. Chattopadhyay and R. R. Sarkar, Journal of Ecological 

Modelling. 163, 45 (2003). 
[13] J. Bascopte and R. V. Sole, J. Theor. Bio. 195, 383 

(1998). 

[14] R. V. Sole and J. Bascopte, Self-Organization in Complex 
Ecosystems, Monographs in Populatin Biology, Princen- 
ton University Press (2006). 



